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In this work, we report on experimental and numerical investigations of the 
attenuation of a small-amplitude signal due to its interaction with high intensity, 
band limited sound whose spectrum consists of one to four discrete peaks. We 
probe the “thermodynamic limit” for different relative frequencies of the spectral 
components. In particular the attenuation of the signal is investigated for both 
equally and unequally spaced spectral components, as well as different phase 
relations among them. The possibility of collective modes is also explored by 
measurements of the phase change in the signal downstream due to the 
presence of discrete noise. 
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I. INTRODUCTION 



The essential process in any system of finite-amplitude noncoherent 
waves is the irreversible redistribution of energy among different frequencies. 
This gives rise to a variety of interesting and potentially useful phenomena, some 
of which have only recently begun to be fully explored. One such effect is the 
absorption of a signal by nonlinear shockless noise. 

The basic nonlinear interaction of acoustic waves is a three-wave 
resonance, so that waves with frequency fi and fa scatter to produce waves with 
frequencies fa = | fi ± fa I . Furthermore, because these waves are nondispersive, 
resonant interactions are collinear (Rudenko and Soluyan, 1977). For random 
waves in two and three dimensions, energy in the noise components 
redistributes irreversibly by collinear resonant interactions. After a time 
determined by the nonlinearities, a closed system of interacting random waves 
can reach thermodynamic equilibrium, and small fluctuations about this state 
relax exponentially with an attenuation coefficient proportional to the noise 
intensity. A small-amplitude wave that is injected into this noise can be regarded 
as a fluctuation of the noise. Whereas the attenuation of this signal is predicted 
(Westervelt, 1976) and observed (Stanton and Beyer, 1978 and 1981) to be an 
exponential for noise in three dimensions, in one dimension the signal is 
predicted (Rudenko and Chirkin, 1975; Rudenko and Soluyan, 1977) to 
attenuate as a Gaussian in the distance from the source. Experimental 
verification of the Gaussian attenuation has been recently reported (Larraza et 
al„ 1996). 

The Gaussian attenuation in one-dimensional nonlinear acoustics brings 
together two physically significant results. First, because this fluctuation does 
not relax exponentially, we infer that a system of interacting, collinear. 
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nondispersive waves is inherently far off equilibrium. This results from both a 
lack of angle randomization and the fact that all spectral components are in 
resonance with each other in this case. A second aspect of a Gaussian 
attenuation is that it breaks translational invariance because the attenuation 
coefficient is now a function of position. The system therefore has memory. 
Specifically, by taking measurements at several distances from the source of a 
signal, one can calculate both the location and strength of the source. This 
remote detection cannot be achieved if the attenuation is exponential. This is 
physically significant because the equations that describe the evolution of 
nonlinear noise and its interaction with a monofrequency signal are 
translationally invariant. 

Because one-dimensional acoustic noise is inherently far off equilibrium, it 
belongs to the general class of systems that exhibit wave turbulence (Zakharov 
et al., 1992). Collective modes, in which energy fluctuations propagate, have 
been predicted in certain far off equilibrium dispersive systems (Larraza and 
Falkovich, 1993), but there are currently no controlled observations. These 
modes are important, for example, as a mechanism of energy transfer. In 
contrast to other wave turbulent systems, nonlinear acoustic noise in one 
dimension is unique in that it lacks a Boltzmann-type description of interacting 
waves. The existence of collective modes in one-dimensional acoustic noise is 
thus unknown. Larraza et al. (1996) have established a first step in this 
possibility: the first observation of the co-^ power law spectrum that is predicted 
(Kuznetsov, 1970) to occur for fully-developed shockless noise. 

Having established some of the properties of broadband noise, a natural 
path of research is to investigate the effect on a small-amplitude signal of band- 
limited noise whose spectrum consists of discrete peaks. In particular, it is 
natural to ask how many peaks would lead to the same absorption of a signal as 
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a continuous band of noise. By reducing the number of peaks while keeping the 
overall intensity constant, we can investigate the interesting fundamental 
problem of how many degrees of freedom are required for the system to behave 
“thermodynamically” (such that the attenuation is Gaussian). 

In this thesis, we report on experimental and numerical investigations of 
the attenuation of a small-amplitude signal due to its interaction with high 
intensity, band-limited sound whose spectrum consists of one to four discrete 
peaks. We probe the “thermodynamic limit” for different relative frequencies of 
the spectral components. In particular, the attenuation of the signal is 
investigated for both equally and unequally spaced spectral cbmpofiehts, as well 
as different phase relations among them. The possibility of collective modes is 
also explored by measurements of the phase change in the signal downstream 
due to the presence of discrete noise. 

A general theoretical background is presented in Ch. II, while Ch. Ill deals 
with the results of numerical investigations in the absence of viscous and thermal 
losses. The results of the experiments are presented in Ch. IV. In Ch. V we 
present the conclusions derived from the research, as well as recommendations 
for further work. 
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II. THEORY 



The propagation of sound in one dimension in a barotropic fluid is 
described by the particularly simple equation (Blackstock, 1972) 



known as the Riemann equation. Equation (2.1) describes the propagation of 
waves in one direction (toward increasing x), were u(x,t) is the fluid particle 



For an ideal gas, the nonlinear coefficient is p = (1+y)/2, where y is the ratio of 
specific heats of the gas. From Eq. (2.1), the instantaneous velocity of 
propagation of a point in the wave with particle velocity u is c+pu. Thus, points 
of the wave having different particle velocity propagate at different speeds. 

A general implicit solution of Eq. (2.1) can be achieved if we note that a 
point of constant particle velocity u of the wave moves with velocity c+pu. The 
particle velocity at (x,t) is thus equal to the particle velocity at (0,x) where the 
retarded time is 



A general implicit solution of Eq. (2.1) with the x = 0 boundary condition u(0,t) = 
f(t) is thus 



(c + pu)ux = 0 , 



( 2 . 1 ) 



velocity and c^ = (dpfdp)o is the equilibrium value of the square speed of sound. 




(2.3) 
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The velocity at the boundary x = 0 thus determines the velocity at x > 0 through 
the retarded time. 

The simplicity of Riemann’s solution allows the determination of exact and 
asymptotic explicit solutions in the preshock region. Examples include Fubini’s 
(1935) pure tone radiation problem, Fenlon’s (1970) generalization to multiple 
frequency sinusoidal source, Kutznetsov’s (1970) asymptotic spectrum for 
broadband noise, and Rudenko and Chirkin’s (1975) Gaussian attenuation of a 
signal due to interactions with normally distributed random noise. 

The main theme in this thesis is to show that such a simple equation 
describes a rich number of fundamental phenofnena. In this chapter, we present 
the analytic solutions to some special cases. Experimentally, the dimensionality 
is controlled by performing measurements in a long traveling wave tube driven at 
frequencies below the first cutoff. Thus, the effects of wall losses become 
important. We show how to account for wall losses and give the modifications to 
the dissipation-free exact solutions. 

A. SHOCK INCEPTION DISTANCE 

At X = 0, suppose particle velocity Ui occurs at time h, and a greater 
particle velocity U 2 occurs at a later time t 2 . Riemann’s equation (2.1) is invalid 
when the second disturbance overtakes the first, which will eventually occur at 
some distance x > 0. This overtaking yields a discontinuity (infinite slope or 
shock) in u. If ti and t 2 differ by a finite amount, a discontinuity will first occur at a 
distance less than x, because the waveform u(x,t) continuously distorts. We are 
thus led to consider velocity disturbances u and u+du and at x = 0 that are an 
infinitesimal time dt apart. If du > 0, the disturbances will coincide at some 
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distance L, after which Riemann’s equation is invalid. If T is the elapsed time, 
then 



[c + pu]T = L , 


(2.4a) 


[c + p(u + du)](T - dt) = L , 


(2.4b) 


corresponding to the first and second disturbances, respectively. 


Multiplying out 


Eq. (2.4b), substituting Eq. (2.4a), and solving for L, yields 

^ ’‘.f. ' j.’ 




, (o + pf)" 
pdf/dt ■ 


(2.5) 



The quantities f and df/dt are the values of the particle velocity and its slope at x 
= 0, respectively. The expression (2.5) gives the distance at which an initial 
positive slope becomes infinite. (A negative slope never becomes infinite.) For 
any solution to Riemann’s equation to be valid at some distance x, it must be 
verified that x < L for all values of f(t). 

For the propagation of a pure tone, the boundary condition at x = 0 may 
be written as 

u(0) = Uq sincot . (2.6) 

The minimum value of the shock length (2.5) corresponds to cot being an integral 
multiple of 2tz. The minimum distance at which shocking occurs for a pure tone 
is thus 
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(2.7) 



PUqCO 



B. SUPPRESSION OF A SIGNAL BY A PUMP 



The simplest problem of interaction of sound with sound is given by the x 
= 0 boundary condition 



where we focus here on the interaction of a strong pump wave of low frequency 
(Op with a weak signal of high frequency cos- The signal has negligible effect on 
the propagation of the pump. However, the pump will modulate the signal and 
generate sidebands at frequencies (0s±<Op, thereby removing energy from the 
signal. An approximation to Fenlon’s (1970) solution for the evolution of the 



The amplitude of the signal depends only on the particle velocity of the pump, 
the frequency of the signal, and the distance. The signal vanishes at the zeros 
of the Bessel function, at which points the energy has been temporarily 
completely pumped into adjacent sidebands. 




( 2 . 8 ) 



signal is 




(2.9) 
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C. ABSORPTION BY NOISE 



For the case of a weak signal of frequency fo in the presence of finite- 
amplitude broadband noise, the amplitude of the signal attenuates as a 
Gaussian exp(-rx^), where 



r = 



2nYfo 2 



-Tins ’ 



( 2 . 10 ) 



where Urms is the rms velocity of the noise (Rudenko and Chirkin, 1975; Rudenko 
and Soluyan, 1977). Broadband noise may be considered by the x = 0 boundary 
condition 



^noise(®>^) ~ B^^COS^COpt + (Pp) , 



n=1 



( 2 . 11 ) 



where B is the peak velocity of each component of the noise, con are densely 
distributed frequencies, and cpn are randomly distributed phases. The Gaussian 
attenuation thus results as a limiting case of multiple pump waves with random 
phases. When a signal is interacting with one pump, there are regions of space 
where restitution occurs (amplitude increases with distance). This effect is lost 
when the interaction occurs with a large number of pumps with random phases. 
Because of the central limit theorem, the noise described by Eq. (2.11) should 
approach a normal distributed noise for a sufficiently large number N. In this 
limit, a theory formulated in terms of the noise (2.12) should accurately yield the 
Gaussian attenuation. 

By squaring Eq. (2.1 1) and averaging over time, we determine the square 
rms particle velocity of the noise. Solving for B yields 
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B = u 



( 2 . 12 ) 




The approach to the Gaussian attenuation can be probed by increasing N while 
keeping Urms constant. 

D. EFFECTS OF WALL LOSSES 



In the absence of wall losses, a monofrequency weak signal interacting 
with noise attenuates as a Gaussian with attenuation coefficient (2.10). The 
amplitude therefore evolves according to 



dx 



pc' 



IxA , 



(2.13) 



where fo is the frequency of the signal and I = pcurms^ is the intensity of the noise. 

Due to thermal and viscous losses, the intensity of the noise does not 
remain constant as a function of position. Specifically, each spectral energy 
component of the noise evolves according to 

If(x) = If (0) exp(-2a^/f X) , (2.14) 

where lf(0) is the spectral shape at the origin and aV? is the linear attenuation 
coefficient (Kinsler et al., 1982). Eq. (2.14) neglects nonlinear losses due to 
energy transfer to other frequencies, and assumes that the dominant 
mechanisms are thermoviscous losses at the wall. In the shockless regime this 
assumption is very reasonable because, due to nonlinearities, fully developed 
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noise is predicted (Kuznetsov, 1970) to have an high frequency tail. Thus, if 
the noise is originally within a band of frequencies Af = f 2 - fi, as is the case in 
the experiment by Larraza et al. (1996), nearly all of the energy remains in this 
band. 

Incorporating thermal and viscous losses at the wall, we generalize Eq. 
(2.13) to 



where I is the total intensity of the noise. Eq. (2.15) is justified by a multiple 
length-scale approach where variations of the amplitude A occur over a shorter 
length scale than that of the noise intensity I, which holds in the experiment of 
Larraza et al. (1996) because the frequency of the signal is substantially greater 
than the upper frequency of the noise. The first term in Eq. (2.15) represents 
losses due to the presence of noise at location x, and the second term 
corresponds to thermal and viscous losses at the walls. Solving Eq. (2.15) for 
the signal amplitude yields 



which was successfully verified experimentally by Larraza et al. 

In the case of interaction with a single pump in the absence of dissipation, 
Eq. (2.9) implies that the amplitude of the signal evolves according to the Bessel 
equation 




(2.15) 




(2.16) 
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(2.17) 



d^A 1 dA 

zr + ; + 

dx X dx 



PUpCO 



\ J 



2 

A 



= 0 , 



with the boundary condition A = Us at x=0. 

Due to thermal and viscous losses, the amplitude of the pump does not 
remain constant as a function of position but varies according to 

Up(x) = Up(0)exp(-a^x) , (2.18) 

where Up(0) is the amplitude of the pump at the origin. In Eq. (2.18) we have 
again ignored losses due to energy transfer to other higher frequencies and 
assumed that the dominant mechanisms are thermoviscous losses at the wall. 

In order to find the effects of wall losses, we can proceed in a fashion similar to 
the case of broadband noise, and assume that except for a pre-factor exp(-aVfs) 
the amplitude of the signal evolves according to Eq. (2.17), with Up given by Eq. 
(2.18). Such an approximation will hold if the length scale of amplitude variations 
of the signal is predominantly dictated by c^/pupCOs, and if the length scale 1/aVfp 
of amplitude variations of the pump is greater than this value. The solution to 
this new equation with the boundary conditions A=Us and dA/dx=0 at x=0 can be 
obtained numerically. In Chapter IV we compare the amplitude predicted by the 
theory with the measured values of a weak signal interacting with a large 
amplitude pump. 
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III. NUMERICAL SIMULATION 



In this chapter, we present the results of numerical integration of a signal 
in the presence of few components of the “noise." We briefly describe the 
method used to integrate Riemann’s equation for sound propagation in one 
dimension in the shockless regime. We specialize to the cases of one, two, 
three, and four noise components with equal and unequal frequency spacing. 

A. FORWARD PROPAGATION METHOD 



As a numerical convenience, we choose p = 1 and c = 1 so that the 
dimensionless Riemann’s equation (2.1) and the expression (2.2) for the 
retarded time become 




d-)| = 0, 



(3.1) 



and 



X = t - , (3.2) 

1 + u 

respectively. As explained in Chapter II, a general implicit solution of Eq. (3.1) 
with the X = 0 boundary condition u(0,t) = f(t) is given by u(x,t) = f(x). 

Eq. (3.2) determines a set of straight lines in the x-x plane. The slope of 
each line is determined by the value of the velocity at the boundary at a given 
time (wavelet). These lines are known as characteristics, where the parameter x 
is the time base for the source signal with t = x at x = 0. 
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The principle of the forward propagation method stems from the fact that 
the characteristics represent the trajectories of motion of the different wavelets in 
the x-T plane. The velocity of propagation of a wavelet is equal to 1+u where u is 
the assigned value of the velocity of the wavelet at the boundary. In the 
propagation of the initial waveform the first incremental distance produces a new 
time associated with the velocity. While at the origin the spacing between time 
steps is commensurate, once the waves propagate a small distance the time 
between adjacent elements is stretched or compressed according to Eq. (3.2), 
which physically corresponds to distortions of the initial waveform. Appendixes A 
and B show a C-language numerical program that implements the forward 
propagation method. 

B. AVERAGING TECHNIQUE 

Our main emphasis is in determining the evolution of the amplitude of a 
monofrequency signal in the presence of noise. As we forward-propagate the 
boundary condition of a time series with the signal plus noise, sidebands of the 
noise appear about the signal’s frequency. Also, the amplitude of the signal is 
expected to attenuate with distance. It is thus necessary to extract the signal 
from the noise at a distance x from the source. We achieve this by Fourier 
analyzing the waveform at the frequency of the signal, which amounts to time 
averaging the in-phase and quadrature components of the waveform. This 
process is referred to as the “discrete Fourier transform” (DFT) because the 
integration is performed with a finite (rather than continuous) time increment and 
is performed over a finite (rather than infinite) time interval. 

The integration in our DFT is accomplished with trapezoidal summation. 
(Simpson’s rule cannot be employed because the time increment is not 
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constant.) Appendixes A and B show the numerical program which includes the 
trapezoidal summation. 

The integration over finite rather than infinite time intervals causes errors 
in the signal amplitude output. For equally spaced noise components, this error 
can be eliminated completely by averaging over an integral number of cycles n of 
the signal such that n = f/Af, where f is the frequency of the signal and Af is the 
spacing between spectral components of the noise (Jang, 1996). However, 
because the time increment is not constant n is only a lower bound. 

We investigated numerically the proper averaging time (measured by the 
cycles of iteration) and its effect on calculated signal amplitude at the fixed 
position x=1 .5 downstream in the system. Figure 3.1 illustrates the result of the 
calculated signal amplitude as a function of number of cycles when the noise is 
composed of four equally spaced spectral components, or pumps, between 1 
and 2 frequency units. The frequency of the signal was chosen to be 50 
frequency units. The relative phase among the four pumps was randomly varied. 
The figure shows that the amplitude of the signal reaches a constant value after 
1,000,000 cycles of the signal. The final output depends on the seed value that 
was used to initialize the relative phases, and thus the signal amplitude output 
depends on the relative phases of the noise components. This result has been 
shown to be true for up to 50 equally spaced spectral components of the noise 
(Denardo et al., 1997). 
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Figure 3.1 Numerical result of the nonlinear interaction of weak signal and four equally spaced 

noise components. The amplitude of the signal is measured at the fixed location x=1 .5 
downstream. The different symbols refer to four different seed values which yield four 
different sets of relative random phases of the noise components. The amplitude of 
the signal depends upon the sets of random phases. 

We next investigated the effect on the signal due to unequally spaced 
noise components. The amplitude of the signal is also calculated at the fixed 
location x=1 .5 downstream. The seed values were used to initialize both the 
relative phases of the noise components, and the unequal frequency spacing. 
The unequal spacing between spectral components of the noise was established 
by adding a random step of roughly 1% to the equally spaced noise components. 
The numerical computer program with the modifications is listed in Appendix A. 
We investigated numerically the proper averaging time in this case as well. A 
number of cycles greater than 1 ,000,000 was needed to get a stable output. 
However, in contrast with the equally spaced case, changing the seed value did 
not change the stable amplitude output, as shown in Fig 3.2. Thus, when the 
noise spectrum is composed of unequally spaced spectral components, the 
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amplitude of the signal amplitude is independent of the relative phases and 
frequencies of the noise components. 




Figure 3.2 Numerical result of the nonlinear interaction of weak signal and four unequally spaced 
noise components. The amplitude of the signal is measured at the fixed location x=1.5 
downstream. Four different seed values where used where each one initialized both 
the relative random phases and the amount of unequal spacing between the noise 
components. The amplitude of the signal is independent of the seed value. 



C. INTERACTION OF A SIGNAL WITH A SMALL NUMBER OF 
PUMPS 



The main conclusion of the previous section is that the attenuation of a 
weak signal due its interaction with discrete pumps that are unequally spaced in 
frequency is independent of the relative phase of the pumps and the value of 
their frequencies. In this section we investigate the spatial dependence of the 
amplitude for one, two, three and four pumps with unequally spaced frequencies. 

As discussed in Sec. 2.2, the suppression of a signal by a pump is 
formulated in terms of the boundary condition (2.8) at x = 0 
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u(0,t) = Up sin(27ifpt) + Us sin(2nfst) 



(3.3) 



where the pump frequency is fp = 1 and the signal frequency is fg = 50. The 
weak signal has amplitude 10"^ and has little effect on the propagation of the 
pump with amplitude 10"^. However, the pump will modulate the signal and 
generate sidebands at frequencies % ± fp, thereby removing energy from the 
signal. In dimensionless variables, Fenlon’s approximate solution (2.9) for the 
evolution of the signal is 



Figure 3.3 shows a comparison between the numerical integration and the theory 
(3.4). The shock inception distance (2.5) for the pump equals 15.92 units. The 
amplitude was calculated by averaging over 10® cycles of the signal. There is 
excellent agreement between theory and numerical integration. 




(3.4) 
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Fig 3.3 Suppression of a signal with frequency 50 by a pump of amplitude 10'^ and frequency 1. 
The signal’s amplitude is normalized to its value of 10"^ at the origin. The points are the 
result of numerical integration, and the solid line corresponds to the theoretical expression 
(3.4). The amplitude of the signal was calculated by averaging over 10® cycles of the 
signal is in the averaging. 



Unlike the case of one pump, there is no closed-form analytic expression 
for the amplitude of a weak signal in the presence of several pumps. Figure 3.4 
shows the evolution of the amplitude of a signal as a function of distance in the 
presence of two pumps, whose frequencies are 1 and 1.25. Note that there is 
suppression of the signal in this case at a dimensionless distance of about 1.1 
(compare with 0.75 for the case of a single pump). Also plotted in Figure 3.4 
are both, the theory Eq. (3.4) and the Gaussian attenuation theory Eq. (2.13) for 
a pump with the same root mean square amplitude and for broadband noise of 
the same energy of the pumps respectively. 
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Fig 3.4 Suppression of a signal with frequency 50 by two pumps. The rms particle velocity of the 
pumps is 10'^ and the frequencies are between 1 and 2. The signal’s amplitude is 
normalized to its value of 10"^ at the origin. The points are the result of numerical 
integration, and the solid lines correspond to the theoretical expression (3.4) and the 
Gaussian attenuation law (2.13). The number of cycles of the signal is 10®. 



Figure 3.5 shows the amplitude of a weak signal with frequency of 50 in 
the presence of three and four pumps with unequally spaced frequencies 
between 1 and 2. We used the values are 1 , 2 ^'^, 2 ^'^ for the frequencies of the 
three pump, and the values 1, 2 ^'^, 2 ^'^, 2 }'^ for the frequencies of the four 
pumps. We also plot the Gaussian attenuation theory for broadband noise with 
the same energy as that of the pumps. Remarkably, the attenuation of the signal 
due to its interaction with four pumps follow closely the Gaussian attenuation due 
to broad band noise. On the other hand, when the signal interacts with three 
pumps, the restitution effect of the signal amplitude is apparent. As we will see 
in Chapter IV, it is also experimentally observed that the amplitude of the signal 
evolves closely as a Gaussian due to its interaction with four pumps, and the 
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departures from Gaussian behavior are more pronounced v\/hen the signal 
interacts with three pumps. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 



distance 



Fig 3.5 Attenuation of a signal with frequency 50 by three (squares) and four (circles) pumps. 

The rms particle velocity of the pumps is 1 0 and the frequencies are 1 , 2 ''*, 2^^ for 
three pumps, and 1 , 2’'*, 2 "^, 2 for four pumps. The signal’s amplitude is normalized to 
its value of 10"^ at the origin. The points are the result of numerical integration, and the 
solid line corresponds to the Gaussian attenuation law (2.13). The number of cycles of 
the signal is 10® in the averaging. 
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IV. EXPERIMENT 



A. APPARATUS 

As stated in the Introduction, the experiment was conducted in a 
waveguide at frequencies below the first cutoff frequency to ensure one- 
dimensional propagation. Air was chosen as the acoustic medium. Other gases 
serve no justifiable advantage, and liquids are undesirable because their typical 
acoustic impedance is not substantially less that that of the confining waveguide 
material. (Such an acoustic impedance mismatch is required if plane-wave 
propagation is to be restricted to the medium that fills a waveguide.) 

Fig. 4.1 shows a schematic diagram of the traveling wave tube built to 
perform the experiment and initially reported by Dorff (1991) . The tube includes 
seven 10-foot (3.05 m) sections of aluminum pipe with wall thickness 1/4 in. (6 
mm) and inner diameter 1 .96 in. (4.98 cm). For this diameter the first (azimuthal) 
cutoff frequency is 4.04 kHz, below which only plane waves can propagate. At 
one end are two compression drivers (J. B. Lansing model 2450H) connected to 
the tube by a “T” adapter, and at the other end is an anechoic termination. 




60 ft 



Figure 4.1 Traveling wave apparatus to observe the absorption of a monofrequency signal by 

nonlinear noise. One driver generates the signal, and the other the noise. The cross 
section of the tube is circular, with inner diameter 2.0 in. Both the signal and noise 
are plane waves. 
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To minimize reflections from the junctions and to reduce possible 
turbulence due to finite-amplitude effects, the sections of the pipe were smoothly 
joined as follows. The surfaces were only approximately circular and concentric, 
so the outer surface of each end was machined to a common diameter, and the 
inner surface was then machined with a slight (1°) taper opening to a common 
diameter at the end. The sections are connected with collars which provide 
alignment, and pairs of flanges which apply pressure between the contacting 
ends. Pulse measurements show that the reflections from the junctions are 
down at least 40 dB (1% in pressure deviation) for frequencies below the cutoff 
frequency of 4.0 kHz (Dorff, 1991). 

Each of the first six sections of the tube has two threaded microphone 
ports with a uniform spacing of 5.00 ft. (1.52 m). An additional port exists 3.5 in. 
(8.9 cm) from the cavity formed by the drivers, and 4.50 ft. (1 .37 m) from the first 
port of the first section. Data from a microphone at the cavity is not useful 
because standing waves prevent an accurate measurement of traveling wave 
amplitudes. Each port accommodates a piezoresistive pressure transducer 
(Endevco model 851 OB) with a 10-32 mounting thread and a 0.15 inch (3.8 mm) 
diameter end that is flush along the inner wall of the tube. We used two such 
matched microphones in the experiment: one remained at the port nearest the 
drivers, and the other was moved from port to port as we gathered data. Each 
port not in use was closed with a screw that was machined to act as a substitute 
for the microphone. 

The anechoic termination spans the length of the last section of the tube. 
The termination consists of 0000-grade steel wool with a tapered linear density 
that increases from zero to a maximum value in a manner according to Burns 
(1971), and is thereafter uniform. The taper occupies 1.0 m of the section; the 
remaining 2.0 m is uniform. The taper was achieved by first piecing together 
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pads of the steel wool to make a large bed, which were then cut according to the 
desired linear density and then rolled into a cylinder. The cylinder was fit into the 
tube by temporarily winding it with string and applying tension. As recorded by 
Dorff (1991), pulse measurements yield reflected waves that are at least 40 dB 
down for frequencies of 500 Hz and greater. 

Tests on sum and difference generation of two signals reveal that the 
nonlinear interactions occur in the medium of interest (air), and are not the result 
of intermodulations of the drivers (Dorff, 1991). This indicates that the apparatus 
is suitable for an experiment to observe the absorption of sound by noise. 

In the experiment, one driver was used to generate discrete “noise” of up 
to four spectral components (referred as pumps), and the other a 
monofrequency signal. The pumps were generated by an HP 8904 
multifunction synthesizer that can provide up to four outputs of different 
frequencies, amplitudes, and relative phases. The power spectral density of the 
pumps was measured with a signal analyzer (HP 3561A). The frequency of the 
signal was typically 3.5 kHz. Due to both the cascading of the noise to higher 
frequencies and the attenuation of the signal, nonnegligible random variations of 
the signal occurred. Hence, to accurately extract the amplitude of the signal 
from the noise, we used a digital lock-in amplifier (SRS 850). The digital lock-in 
amplifier was set with a filter of 24 dB/octave and, for unequally spaced 
frequencies, a time constant of 30 s. Typically, after ten time constants a stable 
amplitude can be recorded. Thus for a 30 s time constant, the averaging time is 
5 min, which for a frequency of 3.5 kHz, corresponds to 10® cycles. Not 
coincidentally, numerically we also required to average over 10® cycles of the 
signal for an stable output. Thus, a digital lock-in amplifier and a digital 
computer agree on the requisite number of cycles for a stable average. 
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B. MEASUREMENTS 



We first established the maximum intensity of the discrete noise within the 
waveguide before shocks were formed. We achieved this by examining a single 
pump generated by the HP 8904. By adjusting the amplitude intensity of the 
pump, and monitoring a microphone output on an oscilloscope, discontinuities in 
the waveform were observed. The waveform was sampled at the final port, just 
before the anechoic termination of the waveguide. An input voltage was 
selected that did not cause the discontinuous waveforms at the final port; thus, 
the noise is ensured to be shockless over the entire length of waveguide. 

For multiple pumps of different frequency, the input voltage of each 
component was determined by taking the maximum input voltage as described 
above, and dividing it by Vn , where N is the number of pumps employed to 
characterize the noise. In this way the total intensity of the noise is constant. 

The importance of maintaining the same noise intensity is to be able to 
accurately compare the effects on the signal due to different number of noise 
components. The various input voltage used in the measurements are listed in 
Table 4.1. 



Number of Components 


1 


2 


3 


4 


Input Drive Voltage 


165 mV 


117 mV 


96 mV 


83 mV 


Noise Spectral Intensity 


480.5 mV 


481.0 mV 


479.8 mV 


479.6 mV 


Noise Spectral Level 


142.45 dB 


142.46 dB 


142.44 dB 


142.43 dB 



Table 4.1 Input drive voltages from the HP8904 for various number of pumps used in 
experimental measurements. 
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The excess attenuation of the signal due to the nonlinear interactions with 
one, two, three, and four pumps was investigated. The frequencies of the 
pumps were selected within the band 0.5 to 1 .0 kHz. As it will become clear 
later, the selection of this band was made because of the need to compare with 
earlier measurements of broadband noise with the same band, conducted by 
Larraza et al. (1996). 

Figure 4.2 shows the attenuation of 3.5 kHz signal by one (fi = 750 Hz), 
two (fi = 633 Hz and fa = 814Hz), three (fi = 633.3 Hz, fa = 762.5 Hz, h = 893.1 
Hz,), and four (fi = 533 Hz, fa = 676 Hz, fa = 780 Hz, and U = 893 Hz) pumps as a 
function of distance. The signal amplitude is the ratio of amplitude in the 
presence of the pumps to that of the amplitude without any pump, thus removing 
effects of wall losses on the signal. 




Figure 4.2 Attenuation of 3.5 kHz signal by 1(open circles), 2 (open triangles), 3 (open 
diamonds), and 4 (open squares) pumps as function of distance. The data 
corresponds to the ratio of the amplitude in the presence of pumps to the amplitude 
without any pump. 
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In agreement with the numerical results of the previous chapter, for 
unequally spaced frequencies of the pumps, we observed that the attenuation of 
the signal is independent of the relative phases of the pumps. However, with 
equally spaced frequencies the attenuation of the 3.5 kHz signal has a strong 
dependence upon the relative phases, also in agreement with the numerical 
results. We verified this statistically with different sets of random phases for both 
three and four pumps. Figures 4.3 and 4.4 shows histograms of the amplitude of 
a 3.5 kHz signal at a distance of 11. 7 m from the source due to interactions with 
three and four pumps, respectively, where the frequencies of the pumps were 
equally spaced. The statistical samples correspond to fifty randomly selected 
sets of relative phases of the pumps. The mean amplitude of the signal at the 
location is 1.0786 mV for three pumps, and 1.266 mV for four pumps. 
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Figure 4.3 Histogram of the amplitude of the signal at 1 1 .7 m from the source for 50 sets of 

relative phases of 3 pumps with equally spaced frequencies: 633 Hz, 763 Hz, and 893 
Hz. The mean value of signal’s amplitude is 1.0786 mV. 
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Figure 4.4 Histogram of the amplitude of the signal at 11 .7 m from the source for 50 sets of 

relative phases of four pumps with equally spaced frequencies: 530 Hz, 675 Hz, 820 
Hz, and 965 Hz. The mean value of signal’s amplitude is 1.266 mV 



We then randomly "tweaked" the equally spaced pumps by a tenth of a 



Hz, causing the spectral components to become unequally spaced. In this case, 



the amplitude of the signal at the same location of 11 .7 m did not change for 



different values of the relative phases. Therefore, the excess attenuation of a 



signal is independent of the variation of the relative phases of the pumps when 



the spectral components are unequally spaced. Moreover, the measured value 



of the amplitude of the signal in the presence of three pumps was 1.07 mV while 



it was 1.13 mV for the case of four pumps. 



Due to wall losses both the signal and the pumps undergo attenuation. A 



quantitative comparison between a theory and experiment requires that the wall 



losses be taken into account. Figure 4.5 show measurements of the amplitude 
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of the signal due to wall losses. The attenuation constant is obtained by a linear 
fit to the data and yields the value 1.256x10'^ Np/(mVHz). This experimentally 
determined value is 7% greater than the theoretical value predicted for our 
apparatus. This difference is mainly a result of the tube not being perfectly rigid. 




Figure 4.5 Attenuation of a 3.5 kHz signal due to wall losses. The attenuation constant a=1.256x 
10'^ Np/(mVHz) was determined from a linear fit to the data. 

Figure 4.6 shows the comparison between the theory, Eqs. (2.17) and 
(2.18), and experiment for the case of a single pump. The differential equation 
(2.17) was numerically solved (Appendix C) with boundary conditions A=1 and 
dA/dx=0, at x=0, and using the experimental parameters fs =3.5 kHz for the 
signal, fp =750 Hz for the pump, Up(0) = 0.64 m/s, and the empirical attenuation 
coefficient a = 0.001256 Np/(mVHz). 
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Figure 4,6 Excess attenuation a 3.5 kHz signal as a function of distance due to the nonlinear 
interaction with one pump. The solid line corresponds to theory, Eqs. (2.17) and 
(2.18), with no adjustable parameters. The dashed line is the theory in the absence of 
wall losses. 



The theoretical curve shows reasonable agreement with experiment. The 
deviations can be understood by noting that the length scale of amplitude 
variations of the pump 1/aVfp « 29 m, while the length scale of amplitude 
variations of the signal due to the pump is c^/(3UpC0s » 7 m or 25% of the length 
scale of the pump’s amplitude variations. Thus, the assumptions for a multiple 
length scale approach may be stretched in this case. 

In the case of four pumps, we can compare the excess attenuation of the 
signal to that of earlier measurements of broadband noise. The broadband noise 
in Larraza et al. (1996) experiment was within a bandwidth of 0.5 to 1.0 kHz and 
had a maximum noise spectral level of 141 dB. The four discrete noise 
components used in our experiment were all within the same bandwidth, were 
not equally spaced, and had a noise spectral level of 143 dB. Figure 4.7 shows 
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the noise spectra for the two cases. In both cases, the spectra is measured at 
the first microphone port, which is 8.9 cm from the drivers. 




Figure 4.7 Comparison of the noise spectrum level (dB re 20 uPa / -v/Plz ) of broadband 

shockless noise levels at the first microphone port, which is 8.9 cm from the drivers, 
to that of the spectrum level of four pumps. The 3.5 kHz signal is also displayed for 
reference. 

Figure 4.8 shows the excess attenuation of the signal for the two cases. It 
should be noted that, due to the slightly higher noise level of the four-component 
case, the signal is expected to exhibit more attenuation. The general 
observation is that four discrete components roughly yield the same nonlinear 
attenuation as that of the broadband noise. 
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Figure 4.8 Excess attenuation of a 3.5 kHz signal due to nonlinear interactions with broadband 
noise (open circles) and with four pumps (open up triangles). The noise band is 
within 0.5 to 1 .0 kHz with an intensity level of 141 dB. The four pumps are within the 
same bandwidth but with frequencies not equally spaced and an intensity of 143 dB. 
The data corresponds to the ratio of the amplitude in the presence of the noise or 
pumps to the amplitude without noise or pumps. The curve corresponds to theory, 
Eq. (2.16), with no adjustable parameters. 



For noise with a discrete number of components, the solution to Eqs. 
(2.15) and (2.16) is 



Ac 



= exp 






'ims 



N ^ 

z 

i=1 V 



AaH 





If we specialize Eq (4.1) to the cases where there are only three and four 
noise components, we obtain the theoretical curves shown in Fig. 4.9. The 
experimental parameters were Urms = 0.6408 m/s for the rms amplitude of a noise 
component, a = 1 .2561 x1 0"^ Np/(mVHz), fg = 3.5 kHz for the frequency of the 
signal. For the 3 pump case, the following frequencies were used fi =633.3 Hz, 
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f 2 = 762.5 Hz, fs = 893.1 Hz, and fi =533 Hz, f 2 = 676 Hz, % = 780 Hz, and U = 
893 Hz for the frequencies of the four unequally spaced pumps. There is a 
surprisingly very good agreement between experiment and the predicted 
Gaussian attenuation for the case of four pumps, but the agreement is not as 
good for three pumps. This indicates that the four unequally spaced pumps may 
be the minimum number of noise components to achieve the “thermodynamic” 
limit of broad band noise. 




Figure 4.9 Excess attenuation of a signal in the presence of three (squares) and four (circles) 
pumps with unequally spaced frequencies. The curves are the theory Eq. (4.1) with 
no adjustable parameters: Dashed curve corresponds to the theory with three 
pumps, and the solid curve is the plot for the theory with four pumps. 



We also conducted a preliminary experiment to examine the phase shift of 
a weak signal due to its interactions with 2, 3, and 4 pumps. The phase was 
measured at two locations, and phase differences were tabulated. The pumps 
had equally spaced frequencies and zero relative phases. The results of these 
measurements are shown in Table 4.2. 
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Number of 


Temp (“C) 


Time 


phase (j> 


phase (j) 


time 


Remarks 


components 






at port 0-B 


at port 4-3 


constant t 





4 


20.5 


1442 


-124.55° 


-101.63° 


1 sec 


without 






14Apr97 








noise 


4 


20.5 


1447 


-124.55° 


-97.54° 


30 sec 


with noise 






14Apr97 










4 


20.5 


1448 


-124.54° 


-98.88° 


1 sec 


without 






14Apr97 








noise 



3 


20.8 


1209 


-124.48° 


-68.39° 


1 sec 


without 






25Apr97 








noise 


3 


20.8 


1214 


-124.46° 


-66.20° 


30 sec 


with noise 






25Apr97 










3 


20.8 


1215 


-124.48° 


-69.02° 


1 sec 


without 






25Apr97 








noise 



2 


20.8 


1247 


-125.48° 


-72.40° 


1 sec 


without 






25Apr97 








noise 


2 


20.8 


1252 


-125.37° 


-70.00° 


30 sec 


with noise 






25Apr97 










2 


20.8 


1253 


-125.47° 


-72.99° 


1 sec 


without 






25Apr97 








noise 



Table 4.2 Measured phases of the 3.5 kHz signal in presence of various numbers of noise 

components. The microphone at port 4-3 is 1 1 m downstream form the microphone at 
the reference port 0-B located 0.87 m from the JBL drivers. 
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The phase recordings were taken from two separate SRS 850 digital lock- 
in amplifiers. The “sync” output of the HP 33120 A function generator used to 
drive the 3.5 kHz signal provided the reference frequency for both lock-in 
amplifiers. There was an observed drift in the phase over time, during what 
appeared to be constant conditions. An analysis of the drift did not appear to be 
related to the temperature variation during this period. To reduce the error due 
to drift, we recorded the pure tone signal phase before and after applying the 
noise components, and averaged the two values. 

When only the signal is driven, the lock-in requires a small time constant, 
while a longer time constant is needed when the pumps are turned on. Taking 
the phase of the signal measured at the downstream port for 4, 3, and 2 pumps 
and subtracting the average phase of the signal in the absence of pumps, 
renders a positive phase shift of 2.715°, 2.505°, and 2.695°, respectively. Thus 
the phase of a weak signal becomes more positive in the presence of high 
intensity sound spectrum. In the presence of pumps, an upward phase shift of 
the signal indicates that the speed of sound increases with increasing intensity of 
the pumps. 
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V. CONCLUSIONS AND FUTURE WORK 



A. CONCLUSIONS 

We have experimentally and numerically investigated the attenuation of a 
monofrequency signal in the presence of a small discrete number of high 
intensity low frequency pumps in one dimension. Analogous to the fundamental 
question in the kinetic theory of gases regarding the minimum number of 
degrees of freedom that can represent the thermodynamic limit, we have probed 
the minimum number of spectral components that can represent noise. Because 
the interaction of a weak signal with broadband noise has been predicted 
(Rudenko and Chirkin, 1975) and experimentally measured (Larraza et al. 1996) 
to yield a Gaussian attenuation of the signal, we have used this attenuation law 
as our probe. 

For equally spaced frequencies of the spectral components of the noise, 
the amplitude of a weak signal due to its interaction with the noise depends on 
the relative phase among the components of the noise. For unequally spaced 
frequencies, the signal’s amplitude is independent of the relative phase. This 
behavior persists even for a larger number of discrete components. On the other 
hand, averaging different values of the signal amplitude corresponding to 
different realizations of the phases of equally spaced frequencies yields the 
same amplitude value as in the case of unequally spaced frequency 
components. 

We observed, both numerically and experimentally, that at least four 
components would yield and attenuation law in close agreement with the 
Gaussian attenuation. This leads to the conclusion that the “thermodynamic” 
limit of nonlinear acoustic noise is nearly represented by four spectral 
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components with unequally spaced frequencies. 

Preliminary observations of the change in phase of a weak signal in the 
presence of two, three, and four noise components reveal that the phase 
becomes more positive in the presence of high intensity sound spectrum. An 
upward phase shift of the signal indicates that the speed of sound increases with 
increasing intensity of the pumps. This constitutes the first observation of such 
phase shift. The significance of this result lies in the possibility of observing 
collective modes, in which energy fluctuations propagate (Larraza and Falkovich, 
1993). Although such modes in far off equilibrium noise have been predicted in 
several systems, there are currently no' controlled observations. These hriodes 
are important, for example, as a mechanism of energy transfer. 

B. FUTURE WORK 

Future investigations should test the statistics of the numerical noise and 
quantify the deviations from normal distribution at the origin. Deviations from the 
normal distribution should also be quantified downstream where strong 
correlations may develop due to nonlinearities. A comparison between equally 
and unequally spaced frequencies should be made, both at the origin and 
downstream. Due to nonlinearities the energy cascade to different frequencies 
may cause a significantly different evolution of the statistics of the noise between 
these two. 

Because a substantially small frequency components of the noise are 
required to yield an approximate Gaussian attenuation of a signal, this may be 
importantin applications where an efficient suppression of high frequency sound 
is desired. 
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Measurements of the phase shift of a weak signal are just the beginning 
of a search for collective modes. The phase shift should be measured as a 
function of the intensity, number of components, relative phases, frequency 
spacing, and also for broadband noise. 
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APPENDIX A. FORWARD TIME NUMERICAL PROGRAM THAT 
DETERMINES THE PROPER NUMBER OF ITERATION CYCLES 
FOR STABLE SIGNAL AMPLITUDE 



/* Numerical integration of Riemann's equation */ 

/* for propagation of sound in one dimension */ 

/* File Name: SIZE.C */ 

/* This program simulates the propagation of a CW in the 
presence of one dimensional noise with a flat 
distribution. The signal's amplitude is calculated as a 
function of distance. The noise is discrete and the 
number of components N can be made to vary from one to 
several. Regardless of the number of components, the 
total energy of the noise is kept constant by dividing 
the peak amplitude of each noise component by the 
square root of N. 

This modification utilizes randomly generated (control by 
varying seed values) numbers to initialize the phase and 
unequal spacing of noise components. In doing so the 
proper number of averages or iteration cycles to yield a 
stable signal amplitude output can be determined. 

LAST UPDATE : 13 May 1997, by Mark Lamczyk */ 
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y/ ★ ★ ★ ★ ★ 



PROCESSOR DIRECTIVES 



* * * * * ! 



#include "math.h" 
#include "stdio.h" 
#include "stdlib.h" 



/***** MACRO DEFINITIONS *****/ 

#define pi 3.141592654 

#define fs 50.0 /* frequency of signal */ 

#define us 0.0001 /* dimensionless peak velocity of one of 

noise components */ 



/***** GLOBAL VARIABLES *****/ 



double X, uO, 
time, 
r[100] , 
rr [100] , 
velocity; 

int N; 



/* distance */ 

/* time */ 

/* random phased */ 

/* random nonperiodic term */ 

/* dimensionless particle velocity */ 



/***** SUB routines *****/ 

double retard (void) , initial_velocity (void) ; 



void main ( ) 

{ /*** begin main ***/ 

int i,-j,k,SIZE, seed; 



double R [100 ] , 
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signal. 



/* extracted signal component 
from noise */ 



old time, new time,/* 



dx, /* 

t,dt, /* 

random, /* 

al, a2, dl, d2; /* 

FILE *fp; /* 



/***** Input Parameters 

printf("\n Enter the number 
scanf("%d", &seed) ; 
printf ( " \nThe seed is %d \n 



retarded time */ 
distance step */ 
time and time step */ 
seed for random number 
generation*/ 

variables used for LOCK-IN */ 

initiates the creation of 
output data file */ 

"sir ★ 'jif '*■ ! 

for the seed: " ) ; 

, seed) ; 



printf("\n Enter the number of Noise components (N>0):"); 

scanf("%d", &N) ; 

printf ("\nN is %d \n", N); 

u0=0 . 01*sqrt ( 1 . 0/N) ; /* dimensionless peak velocity of one 

of noise components/ 

/* thus, the noise energy is the 
input u0^2 */ 



printf ("\n Enter a value for SIZE:"); 
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scanf("%d", &SIZE) ; 

printf ("\nSIZE is %d \n", SIZE) ; 

/***** Variables Initialization *****/ 

dl=d2=t=x=0 . 0 ; 
dt=0 .001; 
dx=0 . 1 ; 



for(i=0; i<100; i++) { 

R[i]=0.0; } 

for(j=0; j<l; j++) { /*** begin j loop ***/ 

srand ( seed) ; 

for(i=0; i<N; i++) { /**begin random generation loop**/ 
r [ i] = (2 . 0*pi*rand ( ) /RAND_MAX) ; /* Random phase of 

noise components */ 

rr [i] = (rand ( ) /RAND_MAX) /lOO . 0; /* Random spacing of 

noise components */ 

printf ("%d %.12f %.12f\n",i, r[i], rr[i]); 

} /*** end random generation loop ***/ 

for()c=l; )c<21; k++) { /*** begin k loop ***/ 

time=0 . 0 ; 

- initial_velocity ( ) ; 
retard ( ) ; 
old time=time; 
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for(i=0; i<SIZE; i++) { /*** begin i loop ***/ 

signal=us*sin (2. O*pi*fs*old_time) ; /* signal 

only */ 

time= (i+1) *dt; 
initial_velocity ( ) ; 
retard ( ) ; 
new_time=time; 

t=i*dt ; 

al=signal*cos (2 . 0*pi*f s*t ) ; 
a2=signal*sin (2 . 0*pi*fs*t ) ; 
dl=dl+al* (new_time-old_time) ; 
d2=d2+a2* (new_time-old_time) ; 
old_time=new_time; 

} /*** end i loop ***/ 

R [k] =R[k] + (2 . 0*fs*sqrt (dl*dl+d2*d2) / (fs*dt*SIZE*us) ) ; 

x=k*dx; 

dl=d2=t=0.0; 



} /*** end k loop ***/ 

seed=seed+l . 0 ; 
x=0 i 0 ; 

} /*** end j loop ***/ 
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x=0 . 0 ; 



/***** Print the result *****/ 

for(i=l; i<21; i++) { /*** begin 2nd i loop ***/ 

x= (i-1) *dx; 

fp=fopen ( "thesisl . dat" , "a" ) ; /* Open the output file 

named as thesisl.dat */ 

fprintf (fp, "% . 12f %d %d %d %.12f\n",x, seed, N, 

SIZE,R[i] /l.O) ; 
fclose (fp) ; 

printf ( "% . 12f %d %d %d %.12f\n",x, seed, N, 

SIZE,R[i] /l.O) ; 

} /*** end 2nd i loop ***/ 

} /*** end main ***/ 



y' ★★★ '3k' sub x'outiriBS 
double initial_velocity ( ) 

{ /*** begin initial_velocity sub routine ***/ 

int i ; 

double vel,var; 
var=0 . 0 ; 

for(i=0; i<N; i++) { /*** begin i loop ***/ 
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/* 

vel=uO*sin (2 . 0*pi* (1.0+ (double) (i) / (double) (N) ) *time+r [i] ) ; 

periodic noise condition */ 

/* NonPeriodic noise condition */ 

vel=u0*sin (2 . 0*pi* (1. 0+ (double) (i)/ (double) (N) +rr [i] ) *time+r 
[i] ) ; 

var=var+vel ; 

} /*** end i loop ***/ 

velocity=var+us*sin (2 . 0*pi*f s*time) ; 
return velocity; 

} /*** gnd initial_velocity sub routine ***/ 



double retard 0 

{ /*** begin retard sub routine ***/ 

time=time-x/ ( 1 . 0+velocity) ; 
return time; 

} /*** end retard sub routine ***/ 
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APPENDIX B. FORWARD TIME NUMERICAL PROGRAM FOR 
THE CASE OF UNEQUALLY SPACED NOISE PUMPS 



/* Numerical integration of Riemann's equation */ 

/* for propagation of sound in one dimension */ 

/* File Name; NONPERIODIC . C */ 

/* This program simulates the propagation of a CW in the 
presence of one dimensional noise with a flat 
distribution. The signal's amplitude is calculated as a 
function of distance. The noise is discrete and the 
number of components N can be made to vary from one to 
several. Regardless of the number of components, the 
total energy of the noise is kept constant by dividing 
the peak amplitude of each noise component by the 
square root of N. 

This modification utilizes randomly generated (control by 
varying seed values ) numbers to initialize the phase and 
unequal spacing of noise components. Fixing the 
frequencies an unequal values was also conducted by 
commenting out the random frequency spacing statement. 

The signal attenuation due to the nonlinear interaction 
with the nonperiodic noise components is determined at 
various downstream locations. The simulation model does 
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not take in account the dissipation to wall losses. 



LAST UPDATE : 28 May 1997 by Mark Lamczyk */ 

/***** PROCESSOR DIRECTIVES *****/ 

#include "math.h" 

#include "stdio.h" 

#include "stdlib.h" 



MACRO DEFINITIONS *****/ 

#define pi 3.141592654 

#define SIZE 1000000 /^numerically determine value that 

yields a stable signal amplitude. 

From THESIS. C */ 

#define fs 50.0 /*frequency of signal */ 

#define us 0.0001 /*dimensionless peak velocity of one of 

noise Components */ 



/***** 



GLOBAL VARIABLES *****/ 



double X, uO, /* 

time, /* 

f[5], /* 

r[100], /* 

rr[100], /* 

velocity; /* 



distance */ 
time */ 

fixed frequencies */ 
random phased */ 
random nonperiodic term */ 
dimensionless particle velocity */ 



int N; 
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SUB ROUTINES 



****** j 



double retard (void) , initial_velocity ( void) ; 

void main ( ) 

{ /*** begin main ***/ 

int i,j,k; 

double R [100] , 

signal, /* extracted signal component 

from noise */ 

old_time, new_time, /* retarded time */ 
dx, /* distance step */ 

t,dt, /* time and time step */ 

seed, /* seed for random number 

generation */ 

al,a2,dl,d2; /* variables used for LOCK-IN */ 

FILE *fp; 

fp=f open ( "out2 . dat " , "w" ) ; /* Open the output file named as 

outl.dat */ 



/***** Input Parameters *****/ 

printf("\n Enter the number of Noise components (N>0):"); 
scanf ( "%d" , &N) ; 
printf("\nN is %d \n", N) ; 
u0=0 . 01*sqrt ( 1 . 0/N) ; 
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seed=0 . 0 ; 



dl=d2=t=x=0 . 0 ; 



dt=0 . 001 ; 
dx=0. 1; 



for(i=0; i<100; i++) { 

R[i]=0.0; } 

srand (seed) ; 

for(i=0; i<N; i++) { /** begin random generation loop */ 
r [i]=(2.0*pi*rand() /RAND_MAX) ; /* Random phase of 

noise components */ 
/* rr [i]=(randO /RAND_MAX) /lOO.O; */ 

/* Random spacing of noise components */ 
} /*** end random generation loop ***/ 

f[0]=1.0; f[l]=1.189207115; 

f [2] =1.25992105; f [3] =1 . 4 14 13562; 

for(k=l; k<17; k++) { /*** begin k loop ***/ 

time=0 . 0; 

initial_velocity ( ) ; 
retard ( ) ; 
old_t ime=time; 

for(i=0; KSIZE; i++) { /*** begin i loop ***/ 
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signal=us*sin (2 . 0*pi*f s*old_time) ; /* signal 

only */ 

time= (i+1) *dt ; 
initial_velocity ( ) ; 
retard ( ) ; 
new_time=time ; 

t=i*dt ; 

al=signal*cos (2.0*pi*fs*t) ; 
a2=signal*sin ( 2 . 0*pi*f s*t ) ; 
dl=dl+al* (new_tirae-old_time) ; 
d2=d2+a2* (new_time-old_time) ; 
o 1 d_t ime=ne w_t ime ; 

} /*** end i loop ***/ 

R[k]=R[k]+ (2. 0*fs*sqrt (dl*dl+d2*d2 ) / (fs*dt*SIZE*us) ) ; 
x=k*dx; dl=d2=t=0.0; 



} /*** end k loop ***/ 

i 

x=0 . 0 ; 

/***** Print the result *****/ 

for(i=l; i<17; i++) { /* Begin second i loop */ 
fprintf (fp,"%f %,12f\n",x,R[i]/1.0); 

x=i*dx; 

} /*** end second i loop ***/ 
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f close ( fp) ; 



} /*** end main ***/ 

^'k'k'k'k'k'k'k'k^'k'k'k'k'k'k g £^OUtin0S ^ 

double initial_velocity ( ) 

{ /*** begin initial_velocity sub routine ***/ 
int i; 

double vel,var; 
var=0 . 0; 



for(i=0; i<N; i++) { /* begin third i loop ***/ 

/* vel=u0*sin (2.0*pi*(1.0+{ ( (double) (i) / (double) (N) ) +rr [i] ) ) 
*time+r[i]); */ /* For Random Freq. Generation */ 

vel=u0*sin ( (2 . 0*pi* (f [i] ) *time) +r [i] ) ; /* Fixed Freqs */ 

var=var+vel ; 

} /*** end third i loop ***/ 

velocity=var+us*sin (2 . 0*pi*f s*time) ; 
return velocity; 

} /*** end initial_velocity sub routine ***/ 



56 



double retard 0 

{ /*** begin retard sub routine ***/ 
time=tiine-x/ ( 1 . 0+velocity ) ; 
return time; 

} z*** end retard sub routine ***/ 
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APPENDIX C. NUMERICAL SOLUTION TO THE MODIFIED 
BESSEL DIFFERENTIAL EQUATION 2.17 



/* Thesis Simulation: Numerical solution to modified */ 
/* bessel differential equation. For the case of Signal */ 
/* attenuation due to the nonlinear interaction with a */ 
/* single noise pump. */ 
/* Naval Postgraduate School, Monterey, California */ 
/* File ‘Name: bessell'c */ 
/* Date last modified: 10 April 1997 by Mark A. Lamczyk */ 



/***** PROCESSOR DIRECTIVES *****/ 

#include "math.h" 

#include "stdio.h" 

#include "stdlib.h" 



/***** MACRO DEFINITIONS *****/ 

#define pi 3.141592654 



/***** GLOBAL VARIABLES *****/ 

float A, B, alpha, x, dx, xmax, amp, aa, ab, ALPHA; 

float fp, fs, beta, c, C, ALPHAstar, Uo, signal; 



main ( ) 

{ /*** begin main ***/ 
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Experimental Constant 



★ ★ ★ T*r ! 



A=7.19; /* boundary condition */ 

B=0.0; /* boundary condition */ 

dx=0.01; xmax=20; 
c=341; beta=1.2; 
fp=750; fs=3500; 

/***** Input Parameters *****/ 

printf ( "Enter a value for alpha:"); 
scant ("%f", Salpha); 
printf ( "\nalpha is %f", alpha); 

printf ("\n Enter a value for input noise intensity ( 

mVolts) : " ) ; 

scant ("%f", ^signal); 

print f (" \nlnput Noise intensity is %f mVolts \n", signal); 

/***** Convert Raw Experimental Constant into 

Theoretical 

Terms *****/ 

Uo= (signal/200.0) *(1.0/61.61)* ( 1 . 0/14 . 7 ) * ( 1 . 0e5/l . 0 ) *(1.0/41 
5.0) ; 

ALPHAstar=-2 . 0*alpha*sqrt ( fp) ; 

ALPHA=alpha*sqrt (fs) ; 

C= (4 . 0* (pi*pi) * (fs*fs) * (beta*beta) *2.0* (Uo*Uo) ) / (c*c*c*c) ; 
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for(x=0.1; x<xmax; x+=dx) { /*** begin x loop ***/ 

/* A & B are the coupled terms of the differential 
equation */ 

B=B-( (1.0/x)*B+(C*exp(ALPHAstar*x) ) *A) *dx; 

A=A+ (B*dx) ; 

aa=A*A; ab=sqrt(aa); /* taking absolute value of A */ 
amp=ab/7.19; /* removing the attenuation due to wall 

losses in order to extract only the 
excess attenuated signal due to the 
presence of the noise */ 
printf("%f \t%f \n" , X, amp) ; 

} /*** end X loop ***/ 

} /*** end main ***/ 
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